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We study the unitary dynamics and the thermalization properties of free-fermion-like Hamilto- 
nians after a sudden quantum quench, extending the results of S. Ziraldo et al, arXiv: 1206.4787 
and PRL (2012). With analytical and numerical arguments, we show that the existence of a sta- 
tionary state and its description with a generalized Gibbs ensemble (GGE) depend crucially on 
the observable considered (local versus extensive) and on the localization properties of the final 
Hamiltonian. We present results on two one-dimensional models, the disordered ID fermionic chain 
with long-range hopping and the disordered Ising/XY spin chain. We analytically prove that, while 
time-averages of one-body operators are perfectly reproduced by GGE (even for finite-size systems, 
if time-integrals are extended beyond revivals), time-averages of many-body operators might show 
clear deviations from the GGE prediction when disorder-induced localization of the eigenstates is 
at play. 
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I. INTRODUCTION 

The concept of ergodicity is at the core of classical sta- 
tistical mechanics: it establishes a connection between 
long-time averages of observables and statistical ensem- 
ble averages- . The extension of the ergodic theorem to 
quantum mechanics was pioneered by von Neumann^ 
in a seminal paper on the unitary dynamics of closed 
quantum systems. The experimental possibility of study- 
ing the nonequilibrium dynamics of "virtually" isolated 
quantum systems - most notably cold atomic species in 
optical latticed^ - has stimulated new interest in this 
issue. The simplest setting for such a nonequilibrium sit- 
uation is that of a sudden quench of the Hamiltonian, 
Ho ^ H, governing the dynamics of the isolated quan- 
tum system: the after-quench (t > 0) unitary evolution 

is simply |*(t)) = e^^^*\'^'o) but the initial state {"^q) 
— for instance, the ground state of Hq — is generally a 
complicated superposition of the eigenstates of H. Since 
the energy is conserved, {^{t)\H\'^{t)) = (*ol-H^I*o), it 
is reasonable to expect that a "generic" ergodic evo- 
lution in the Hilbert space will lead to time-averages 
which are reproduced by the micro-canonical ensemble: 
this is wh at v on Neumann discussed for "macroscopic" 
observable^^m^ and is generally expected to occvc^^, in- 
dependently of the initial state |^'o). 

Classical physics teaches us also when to expect viola- 
tions of ergodicity: on one extreme, for systems that are 
integrabl^ or close-enough to being integrabl^^Smi on 
the other, for systems with a glassy dynamics, be it due 
to interactions providing dynamical constraintJ^^]^ or to 
genuine disordeJ^^. Quantum mechanically, dangers to 
ergodicity come from very similar sources: integrability, 
interactions, and disorder. Integrability implies the ex- 
istence of many constants of motion J^, and this clearly 
restricts the ergodic exploration of the micro-canonical 



energy shell, leading to what one might call a "break- 
down of thermalization'^4ll5J_ often results in a kind 
of "generalized thermalization" described by a statisti- 
cal ensemble which maximizes entropy in the presence 
of the constraints (/^), an ensemble introduced long-ago 
by JayneJi^ and known as generalized Gibbs ensemble 
(GGE) IHHS] 



Pgge 



Tr e 



(1) 



where are Lagrange multipliers which constrain the 



mean value of each to its t 
Tr ' 



value: 



(*o|/mI*o> 



On the theory side, the approach to equilibrium has 
been carefully investigated for one-dimensional Bose- 
Hubbard models describing quench experiments in op- 
tical lattices and super-lattice^^HtSIl^ Starting from 
non- homogeneous initial states (density waves), such 
as l^'o) = |1, 0, 1, 0, . . . , 1, 0), and evolving the system 
with a Bose-Hubbard H = + H.c.) -f 

{U/2) J2i ni{ni — 1) both at the integrable points U = 
(free bosons) and U — oo (free hard-core bosons/spinless 
fermions), where analytic solutions are possibl^^Il^ and 
at general (nonintegrable) values of U (through time- 
dependent DMRCP^^^, the physical picture emerging 
has lead to the so-called local relaxation conjectur^^. 
Although the system is, sensu stricto, in a pure state 
p{t) = \^{t)){'^{t)\, when measured-upon locally in a fi- 
nite region S, the resulting (mixed) reduced density ma- 
trix ps{t) — ^^L\sP{t) relaxes towards a stationary Gibbs 
state of maximum entropy compatible with all the con- 
stants of motion. This relaxation is strongly tight to the 
"locality" of the observable O one measures, and results 
from "information transfer" carried by the excitations 
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along the systenPHMlS]^ which eventually "thermalizes" 
any finite region S, the rest of the system acting as an 
effective bath for S. Remarkably, such an approach to 
equilibrium does not require time-average^^Il, Experi- 
mentally, a fast dynamical relaxation was recently ob- 
servecP^ in a system of cold atoms, where the long-time 
stationary state was compatible with GGE. 

A far more dangerous type of "ergodicity crisis" de- 
rives, apparently, from dynamical constraints imposed by 
interactions^^: when quenching a Bose-Hubbard model 
starting, for instance, from l^fo) = |2, 0, 2, 0, . . . , 2, 0) — 
i.e, with an average density of bosons per site n — 1 
where, in equilibrium, a Mott phase transition occurs for 
U > Uc ^ 3.5 J — the ensuing dynamics leads to a fast 
relaxation/thermalization for quenches at small U, while 
the relaxation is extremely slow (and the more so, the 
more the size of the system increases) and the dynamics 
appears effectively freezed for U > [/^f^29l_ Concerning 
ergodicity-breaking due to genuine disorder, both Ander- 
son localization, at the single-particle level'^", and many- 
body localization, in the presence of interactions^^, are 
well known examples of disorder-induced phenomena oc- 
curring in equilibrium physics. Quantum quenches in the 
presence of disorder and interactions have also been stud- 
ied, in the framework of many-body localizatiorP^HlIl^ 
the physical picture is far from being fully understood. 

While breaking of translational invariance and disorder 
in the initial state |5'o) have apparently little or no effect 
on the ensuing relaxation towards a stationary state if 
th e after-q uench Hamiltonian H is translationally invari- 
anlPSIlilllHill^ disorder in the final Hamiltonian H seems 
to play a quite crucial role: this has been seen in numer- 
ical studies of nonintegrable fermion models with disor- 
dered long-range hopping and nearest-neighbor interac- 
tiorPS', and has been recently corroborated through an- 
alytical arguments on integrable quantum spin-chains of 
the Ising/XY clasPl. 

Our present study extends that of Rcf. 38 , where it was 
shown that, in presence of disorder in H, the long-time 
after-quench dynamics does not relax towards a station- 
ary state, and time-fluctuations generally persist in the 
expectation values of local operators, even in the ther- 
modynamic limit; this is, essentially, due to the pres- 
ence of a pure-point spectrum of H associated to localized 
wave-functions, as opposed to the smooth continuum of 
a system with extended states. One can view this per- 
sistence of time-fluctuations in local measurements as a 
result of the inability of the system to carry informa- 
tion arouncP^, due to localization. Time-averages are, 
therefore, mandatory in comparing dynamical quanti- 
ties to statistical ensemble averages. But here a fur- 
ther surprise emerges, which will be the main focus of 
this study: while time-averages of one-body local ob- 
servables are perfectly well reproduced by the relevant 
statistical ensemble — the GGE, being the considered 
models, essentially, free-fermion ones — time-averages of 
many-body local operators show clear deviations from 
the corresponding GGE prediction, in agreement with 



recent numerical findingP^. 

We will exemplify these ideas on a model of dis- 
ordered one-dimensional fermions with long-range hop- 
pings, showing a transition between extended and lo- 
calized eigenstates, and on disordered Ising/XY models. 
We will start. Sec. by stating in a more precise way 
the problem we want to analyze. Next, Sec. |III[ we de- 
scribe the models we have investigated. In Sec. |IV| we 
will analyze the essential reason why GGE works per- 
fectly well in predicting infinite-time averages of one-body 
operators, even for finite systems (as long as one inte- 
grates over revivals). In Sec. |v] we will discuss why this 
is not generally the case for many-body operators, unless 
time-fluctuations of one-body Green's functions vanish 
for large t. Next, in Sees. | VI| and | VII| we will present our 
results for the two models we have considered. Sec. IVIIII 
contains a final discussion, and our conclusions and per- 
spectives. 

II. STATEMENT OF THE PROBLEM 

Let us start by defining the problem. Consider a stan- 
dard quantum quench: an initial state |^'o), ground state 
of some Hamiltonian Hq, evolves under a different time- 
independent Hamiltonian H. Given an observable O, its 
quantum average can be separated in two terms (we take 

0{t) = (vl/o I e'^* 06-^^*1*0) = (6)...„. + 50{t) , (2) 

where the time-independent part (O)time results from an 
infinite-time average of 0{t) 

{OU^^ = lim i rdt 0{t) , (3) 

while SO{t) represents a fluctuating part (with a vanish- 
ing time-average). An alternative standard decomposi- 
tion of 0{t) proceeds by introducing the (many-body) 
eigenstates |a) of H, with energy Ea, and inserting reso- 
lutions of the identity in Eq. Defining Cq, = (a|^o) 
and Oa'a = (a'|0|a), one gets: 

Oit) = ^|apO„,+ ^ e^(^'''-^°)*C:,0„,„C„ 

/ + 00 
dne-'^^'Fom, 
-oo 

where the first term, (O)o, is the so-called diag- 
onal averag^^mH^ while the second, time-dependent, 
part has been recast as a Fourier transform of a 
weighted joint (many-body) density of states Fo{^) = 
j:o.'^o.C*^'Oc.'M{n- E^+E^,). With the quite 
usual assumption of no energy degeneracy, E^' ^ Ea if 
a' a, one can quickly prove that the dia gonal average 
indeed coincides with the long-time averag^^iaifil; 

(0)d = (0)._ = lim ^ fdt 0{t) , (4) 
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while time-fluctuations are given by the Fourier trans- 
form of Foiny. 

/+00 
dne-'^'Fo{n) . (5) 
-co 

As discussed in Ref. UHl the behavior of the fluctuating 
part SO(t), relaxing and decaying to or remaining fi- 
nite (with persistent oscillations) for i — oo, is strongly 
tied to the "smoothness" of Fo{^) in the thermodynamic 
limit (for finite systems, Fo{^) is always a series of dis- 
crete Dirac's deltas, hence SO{t) will never go to zero for 
t — > oo, and revivals will appear). Indeed, SO(t) will de- 
cay to zero for large t if Fo{^) is smooth enough, due to 
the destructive interference induced in the Jl-integral by 
the strongly oscillating phase e~'^* (Riemann-Lebesgue 
lemma); on the contrary, disorder and an important pure- 
point spectrum part, i.e., delta functions associated to lo- 
calized eigenstates which do not merge smoothly into a 
continuum, will lead to persistent time-fluctuations 60{t) 
for local operator&3Si. 

In the rest of the paper we will concentrate on two 
"solvable" models with disorder which will be presented 
in Sec. IIIII disordered one-dimensional fermions with 
long-range hoppings (which can show either power-law 
localized states or extended ones, depending on a pa- 
rameter controlling the long-range hopping variance), 
and disordered Ising/XY model. Both, being quadratic 
fermionic models, can be quite effectively "numerically 
diagonalized" , i.e., one can flnd the one-body spectrum 

and the corresponding quasi-particle creation operator 

7^ to express H = ^iiljilfj.- This, in turn, allows us to 
calculate the fermionic local one-body Green's functions: 

Gnnit) ^ (*(i)l4c,J*(i)) 

F,,,,{t) ^ i^mlcll^it)) , (6) 

where c| creates a fermion at site j. We will show that, 
in view of Wick's theorem, a key issue in understanding 
the validity of GGE for general many-body observables 
in such free-fermion-like Hamiltonians has to do with 
the nature of the long-time fluctuations of the one-body 
Green's functions. That might seem a simple matter to 
explore, but unfortunately, in the general disordered case, 
the one-body Green's functions are essentially impossible 
to obtain analytically, and their numerical study is often 
elusive: you can only study a finite system-size L for a 
finite time t, and whether the fiuctuations will eventu- 
ally vanish or not for L — > oo (first) and i — oo (after, 
otherwise you always get revivals) is often hard to tell. 
It turns out that a useful tool to distinguish the presence 
or absence of time fluctuations for large t is given by the 
time-averaged fluctuations: 

4^ hm ^ rmt)f , (7) 

which is zero when the fluctuations vanish for large i, and 
is flnite otherwise. We will show how to make progress 



analytically, for this quantity, to pin-point the behavior 
of the large-t fluctuations of one-body Green's functions. 
This analytical progress will require an assumption of ab- 
sence of gap degeneracies, i.e., e^^ — e^^ = e^^ — e^^ only 
when /ii — fj,2 and /is = /i4 or /ii — /is and fJ,2 = ^J'4, re- 
stricted however to the one-body spectrum, an assumption 
that can be argued to be a reasonable one, especially in 
presence of disorder. Its general validity for many-body 
energy eigenvalues, on the contrary, is rather tricky, see 
discussion in Sec. IVIIII 

In order to get smoother results we will perform av- 
erages of Sq over different realizations of disorder. We 
stress, however, that the disorder average will be al- 
ways performed after the computation of 6q for each 
realization: performing the disorder average before tak- 
ing the squared time-integral, i.e., on SO{t), would ef- 
fectively kill the time- fluctuations SO{t), by a kind of 
"self-averaging"'SI, 



III. MODELS 

In this work we concentrate on two "solvable" mod- 
els possessing a simple fermionic description. The flrst 
model describes spinless fermions hopping on a chairP^l. 

Hhop = -^nh (c]iCj2 + H.c.) , (8) 

where (cj) creates (destroys) a fermion at site j and 
Jjij2 is a (real) hopping integral between sites ji and 
j2. We will in general take the Jjua's to be random and 
long-ranged, with a Gaussian distribution of zero mean, 
(Jjjjj) = 0, and variance given by: 



Here a is a real positive parameter setting how fast the 
hoppings' variance decays with distance. The peculiar- 
ity of this long-range-hopping model is that, regardless 
of the value of (3 (which hereafter is fixed to 1), it has 
an Anderson transition from (metallic) extended eigen- 
states, for a < 1, to (insulating) power-law localized 
eigenstates for a > Physically, this is due to the 

fact that, for small a, long-range hoppings are capable of 
overcoming the localization due to disorder. The clean 
nearest-neighbor hopping model is recovered by taking 
"/?ij2 = ~^i2,ji±ij where Sij is the Kronecker's delta: we 
will always use this choice for the initial Hamiltonian Hq, 
with the corresponding ground state I'i'o) being the filled 
Fermi sea. (The reason behind this simple choice for Hq is 
that the long-time fluctuation properties do not depend, 
qualitatively, from the initial Hamiltonian being ordered 
or not, see Ref. f55|) . Being quadratic in the fermion op- 
erators, Hfiop can be diagonalized for any chain of size L 
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in terms of new fermionic operators 



(10) 



where Uj^ are the wave-functions of the eigenmodes of 
energy e^: Hhop = Z)^ <^mcJiC^- The energies e^, and the 
associated wave-functions are obtained, for any given 
reahzation of the hoppings Jjijj in a chain of size L with 
open boundary conditions, by numericahy diagonahzing 
the L X L one-body hopping matrix. 

The second Hamiltonian we considered describes a dis- 
ordered Ising/XY chain in a transverse fielcP^ 



H 



XY 



L 

E 



(J-ap;Vi + ^-J-]Vi)-EVf, (11) 
i=i j=i 



where L is the size of the chain, ct^ (fi = x, y, z) 
are spin-1/2 Pauh matrices for the j site, and periodic 
boundary conditions are assumed, =0"^. JJ, Jj 

and hj are real and describe, respectively, the nearest- 
neighbor spin couplings and the transverse magnetic 
field. A quadratic fermionic Hamiltonian is obtained 
here by applying a Jordan- Wigner transformation"*® c/ = 

exp (^iTrX^jCl ^t^J^ ' where af = (ctJ ± ia^)/2: 



while, for the final disordered Hamiltonian, we take Jj = 
1 + er]j ,7 = 1 (Ising case) , and hj = h + e£,j , where e sets 
the disorder strength and rij , are uncorrelated uniform 
random numbers in [— 1, 1[. 

The aim of our analysis is to study the behavior of 
different disorder realizations for very large L (ideally, in 
the thermodynamic limit). Practically, the largest size 
we will consider is L = 2048. Given a disorder realization 
for L = 2048, we will generate corresponding realizations 
at smaller L by cutting away the same amount of sites 
from the two edges. In this way, we obtain a smoother 
behavior for all quantities versus L. For every quantity x 
considered, we have checked its probability distribution 
P{x) for different realizations of disorder. Sometimes 
P{x) deviates strongly from a Gaussian distribution, and 
is nearly (although not precisely) log-normal, i.e. it is 
log(a;) which is approximately Gaussian-distributed. In 
such a situation, we will calculate and plot the median 
(i.e., the geometric mean) [x]^ = exp ( [log a;]^^^) and the 
geometric standard deviation exp(fT[log x]), rather than 
the usual (arithmetic) mean [x]^^, and its standard devi- 
ation a[x] (in the plots, the error bars will then go from 
Wavexp(-cr[logx]) to [x]^^ exp(cr [log x])) . 

For both the models considered above, after diagonal- 
ization, the Hamiltonian is expressed as: 



H 



(14) 



HxY= - ^-^(4 

L 

- E ^^(246,-1), 



Cj+i + 744+1 + 



(12) 



where Jj = Jj + Jj, and the anisotropy parameter 7 is 
such that JJ = Jj{l + j)/2 and Jj = Jj{l--f)/2. In the 
relevant sub-sector with an even number of fermions, we 
have to apply anti-periodic boundary conditions c^+i = 

.m 



where is the positive excitation energy of the state 
7^|0) and Eq is the energy of the state |0) annihilated by 

all the 7^. Notice that, in diagonalizing Hhop, Eq- (|8|, 
one generally obtains some negative e^: in such a case, it 
is enough to perform a particle-hole transformation 7^ — 
cji to change the sign of e^. Physically, that implies that 
all negative energies are occupied in the ground state |0), 
and the resulting excitations describe particles or holes. 
In the following two sections we will use expression 



( 14 1 for H: all the observations made are valid for both 



c^L^ At variance with the fermion- hopping case, thFmodels we have just described. Notice that the num 



Eq. ([8|, there are now, for 7^0, BCS terms 44-i-i which 
create (and destroy) pairs of fermions. Using a Bo goli- 
ubov rotation we therefore define new 

fermionfpaail 



ber operators 7^7^ commute with H in Eq. ( 14) and are 



L 



(13) 



which diagonalize the Hamiltonian Hxv ~ 
I]^=i (7^7^ - 1/2)- The (positive) eigenvalues 
e^/2 and the associated eigenfunctions {uj^^Vj^) are 
obtained, once again, by diagonalizing a 2L x 2i 
one-body matrijP^SSl. Notice also the strict particle- 
hole symmetry present even in the general disordered 
cas^^^HSl for every positive eigenvalue e^/2 > 0, with 
associated (uj^,Wj>), there is a negative eigenvalue 
— e^/2 associated to {v*^,u*^). 

When considering quenches for this Hamiltonian, we 
always start from a clean Hq, with Jj = 1 and hj = ho, 



therefore obvious constants of motion in the GGE aver- 
ages, 4 = 7^%- 



IV. WHY GGE WORKS FOR ONE-BODY 
OBSERVABLES 

In this section we are going to show that, for a free 
fermion Hamiltonian of the form given in Eq. (14), the 



long-time average (and the diagonal average) of any one- 
body operator coincide with the corresponding GGE av- 
erage, for any system size L, and any possible quench: 



{Ol-body) tin-io — {Ol-body)-D — {Ol-body)c 



(15) 



with some important remarks to be made when there are 
degeneracies in the one-body spectrum (see below). As 
we shall see, this can be traced back to the constraints 
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that the GGE sets through the constants of motion 
We stress that, remarkably, this equahty holds even for 
a finite size chain L, while, usually, statistical ensembles 
need a thermodynamic limit. Examples of one-body op- 



erators are, in real space, cj^ or , the local density 
c^jCj, the density fi = ^-fij/L, and, in momen- 



tum space, c|,Cfc (where cl ^J^i 



k - L.j^"^'c]/^/L), etc. More 
generally, a one-body fermionic operator can always be 
written, neglecting irrelevant constants and rewriting the 
Cj's in terms of the 7^ (inverting Eqs. (10) or (13)), as: 



Concerning the equality {Oi-body)ti^.o = {Oi-body)o, 
we should pay attention to the cases in which H has de- 
generate single-particle eigenvalues, e^^ = e^^a^^^i, (for 
instance, when H is disorder- free) . In these cases, the 
time-average of Oi^bodyit) suppresses all the oscillatory 
factors e^''-'^''2+^''i)* occurring in the B— and D— terms 
of Eq. (16), but all the factors e^^'^'^^-'^i^i)'^ correspond- 



fJ.lfJ.2 

(16) 
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where A, B, D are L x L matrices. Let us start showing 
that {di-body)u = (Oi-body) GGE- If I a) = 7^17^2 ' ' ' |0) 
denotes a general many-body eigenstate of H, then 
clearly only the diagonal elements of A enter in the di- 
agonal matrix element: 



{a\d 



ing to degenerate eigenvalues appearing in the A— terms 
survive. Therefore: 

(Oi_body)ti„.o = Y ^MlM2 (*0 17^1 7m2 1*0) 

(/'l,A'2) 

= {dl-body)n+ Y ^MiM2(*0|7ii7/X2l*0) , (20) 

where we have singled-out the diagonal elements, and the 
second sum runs over all the degenerate pairs {fj,i , /i2 7^ 
fj,i) such that — e^^. Because of degeneracies, how- 
ever, there is more freedom in the choice of the fermionic 
operators 7^: we can always perform a unitary rota- 
tion in each degenerate subspace in such a way that 

vt 



1 — body 



c^) = y A^^(a|7t7^|a) , (17) (*o|7ii7p2 l*o) = for m ^ ^2- With such a choice 
^ of the 'yVs, the extra terms in Ec^. (20 1 vanish, and we 



recover the initial statement in Eq. (151, i.e., for any size 



where (ctlT^Jfila) = 0, 1 is the occupation number 
of eigenmode /i in the eigenstate \a). In terms of 
{a\Oi-body\o:) , the diagonal average of Oi-body is readily 
expressed as: 

{dl-body)n = ^^|CapA^^(a|7i7Ml") 

= II^MMlI(*o|a)(a|7i7Mla)(a|*o) 

fj, a 

= E^MM(*o|7i%l*o) , (18) 



where we have added and extra sum over a', using 
that {a\jl^fj,\a') = da,a'{a\jljf,\a), and then recog- 
nized two resolutions of the identity No- 
tice, therefore, that the initial state enters only through 
(^'o|7^7^|5'o), i-e., exactly the constants of motion which 
are constrained and reproduced by the GGE averages: 
(7i7/i)GGE = (*o|7i7;i|*o)- The conclusion is therefore 
simple, as the GGE average of Oi-body is: 

(Ol_f,odl;)GGE = A^fj^ {lulti) GGE 

= 5^A^^(vI.o|7i%|*o) , (19) 
where we used that, by construction of the GGE, 

(7ii7M2#A^l)GGE = {iUiIJgGB = (7mi7m2)gGE = 0. 



and any quench, the long-time average of any one-body 
operator is equal to the GGE one. We stress the fact 
that, for any finite system, Oi-body(t) will have recur- 
rent fluctuations SOi-bodyit), the so-called returns or re- 
vivals, due to the discreteness of the finite-size spectrum: 
nevertheless, integrating over all times (across revivals) 
is guaranteed to reproduce the GGE average: 



lim — 

T->oo T 



dt Oi-body{t) — {Oi-body) GGE 



This statement, however, does not imply that the fluc- 
tuating part SOi^body(t) decreases to for t — oo, as 
indeed evident from the presence of finite-size revivals. 
As we shall see, SOi-body{t) might indeed persist for all 
times even in the thermodynamic limit, when H is disor- 
dered: this in turn implies that limt_>.oo Oi-body{t) might 
not exist in some cases, preventing a straightforward ap- 
plication of Wick's theorem to extend the equalities of 



averages in Eq. ( 15 ) to many-body operators. 



V. GGE FOR MANY-BODY OBSERVABLES 

In the previous section we have shown that the GGE 
average of a one-body operator coincides exactly with its 
long-time average. Here we will see that, for a general 
many-body observable O, the situation is more compli- 
cated, and GGE can be proven to predict correctly long- 
time averages under two additional requirements: 1) O is 
a finite sum of powers of some fermionic operators, and 
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2) the time-fluctuations of the one-body Green's func- 
tions associated to such fermionic operators are vanish- 
ing. Whenever either of the two conditions is not real- 
ized, (O) GGE is not guaranteed to coincide with (O)tin,o: 
we will discuss indeed cases of a definite disagreement 
between the two averages because of a violation of con- 
dition 2) above, i.e., the persistence of time-fluctuations 
of one-body Green's functions. 

The key to the story is Wick's theorerri^, which 
clearly applies to the free-fermion Hamiltonians we are 
discussingSB. Any many-body observable O can be ex- 
pressed as a linear combination of powers of the orig- 
inal real-space fermions Cj and cj. Since |4'o) is a 

BCS-Slater determinant, we can expand {t)\0\^ (t)) , 
using Wick's theorem, as a sum of products of one- 
body Green's functions Gj^j^{t) = {'^{t)\c'-^CjJ^{t)) and 

Fj^j^{t) = (*(i)|c]^c]j^'(t)). To make things more clear, 
let us consider, for instance, the density-density correla- 



tions 



'Jl "-32 



with ji ^ j2j a two-body operator 



(21) 



whose Wick's expansion reads: 

This expansion clearly involves a finite number of terms 
(condition 1). Now suppose (condition 2) that the 
time-fluctuations of the Green's functions vanish for 
large t, hence the limits limt_j.oo Gj^^j^ (t) = Gj-^j^ (oo) 
and limt^oo Fjij.it) ~ Pjihi'^) exist. From the 
analysis of the previous section, it is obvious that 
such limits must coincide with the corresponding GGE 
averages: Gj^j^{(X)) = (c]iCj2>GGE 



_ ^GGE 



It follows there- 



fore that limi_j.oo Pji j2 (^) exists (i.e., its fluctuating part 
^Pjij2 (0 vanishes for large t) and is given by: 

lim n (i^ — r^GGE^GGE |^GGE|2 , I pGGE|2 

~ (PjiJ2)gge , (22) 

where the final step uses the fact that Wick's the- 
orem also applies to GGE averages of free-fermion 
Hamiltonians^. Notice that, since long-time fluctuations 
of Pj^j2 (t) vanish, the infinite-time limit dominates the 
time average, and this implies: 



(Pjij 



lim Pj.j^it) = (p-i.-i^) 



J2/GGE 



For a similar discussion when H is a translationally in- 
variant disorder-free Ising-XX chain, see Ref. 2L A sim- 
ilar proof works quite generally for all observables O 
provided the two stipulated conditions are satisfied. A 
few important remarks are in order: i) The existence 
of time-limits of Green's functions, as opposed to time- 
averages, are crucial in applying Wick's theorem, be- 
cause it is generally false that the "time-average of a 
sum of products" coincides with the "sum of products 
of time- averages"; ii) For definiteness, we have chosen. 



above, the real-space fermions c] to expand O, but sim- 
ilar arguments can be make in any one-body fermionic 
basis for instance a momentum space basis. No- 
tice, in this respect, that O might involve an infinite ex- 
pansion in terms of the c|'s and a finite one in terms 
of the //'s (condition 1). In this case, if the time- 
fluctuations of Glmit) = {'^{t)\flfm\'i'{t)) and J'hnit) = 

(^(i)l///ml^(^)) vanish (condition 2) one can still con- 
clude that (O)timo = (O)gge; Hi) Whenever the time- 
fluctuations of the one-body Green's functions do not 
vanish for large i, and/or the expansion of the operator O 
involves an inflnite number of Wick's contractions there 
is no guarantee that GGE will not reproduce long-time 
averages: we simply cannot prove it by using Wick's the- 
orem. Nevertheless, we will discuss later, see Sec. |VIB 



explicit cases where the persistence of one-body time- 
fluctuations, due to disorder and to the presence of lo- 
calized eigenstates, indeed leads to a definite discrepancy 
between (0)ti„e and (O)gge- 



VI. ID SPINLESS FERMIONS WITH 
LONG-RANGE HOPPING: RESULTS 

Let us consider quantum quenches with a final Hamil- 
tonian given by Hhop, Eq. ([s]), describing disordered 
one-dimensional spinless fermions with long-range hop- 
ping. We will consider observables O in real space and 
in momentum space, and hence we will need to ascertain 
the time-dependence of both real-space and momentum- 
space Green's functions. To unify the treatment of 
both cases, we will consider a general fermionic operator 
/„ = J2f_i UnfiCf^ obtained by applying a unitary transfor- 
mation u (of matrix elements Un/i, with u^u the identity) 
to the c^'s which diagonalize Hhop- For the original real- 
space fermions Cj, Uj^^ is the real-space wave- function of 
the /i-eigenstate, while for the momentum space f erm ions 



IV 



and 



Ck, Ukfj, — 6 '^''^Ujfj,/y/L. As discussed in Sees. 

|vj if the Green's functions associated to the /„'s have 
vanishing long-time fluctuations, then also the long-time 
fiuctuations of O disappear and (O)timo = (O)gge- We 
will show the crucial role played by the localization of 
the eigenfunctions, which we will characterize through 
the standard Inverse Participation Ratio (IPR)P3. When 
quenching to a final Hamiltonian with a in the localized 
phase (a > 1), the real-space Green's functions Gj^j^^t) 
will be shown to have persistent time-fluctuations, i.e. 
(5q. . > in the thermodynamic limit L oo, while 
quenching with a in the extended phase {a < 1) leads 
to vanishing time-fluctuations, 53^ — ^ 0. In both 
cases, however, the eigenfunctions appear to be extended 
when analyzed in momentum space, which results in mo- 
mentum space Green's functions with vanishing time- 
fluctuations, Sq^ ^ — 0. We will then explicitly discuss, 
see Sec. |VIBl the discrepancy between time-averages and 
GGE averages for real-space many-body operators, like 
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density-density spatial correlations, when quenching to 
the localized phase. 



A. One-body Green's functions fluctuations 



As discussed in Sec. IV GGE predicts correctly the 
infinite-time average of any one-body operator, and 
hence, in particular, of the one-body Green's functions. 
Figure [T] exemplifies this, showing the time evolution of 



o 
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Extended phase (q=0..'3) - L = 12& 



- time evolution 
GGE average 
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FIG. 1: Time evolution of Gjj{t), with j — L/2 (solid line), 
and its running-time average t^^ Jq dt'Gjj{t') (dashed line), 
for two different values of a and two different sizes. The 
horizontal dash-dotted line is the GGE average of hj. The 
data are obtained using a single realization, with a = 0.5 for 
the extended case, and a = 2 for the localized one (which 
shows no size effect). Disorder realizations for the smaller L 
chain are obtained, here and in the following, by removing 
sites from the right and left edges of the larger chain. 



Gjj{t), the expectation value of the fermion density at 
site j, for two sizes and two values of a, one in the ex- 
tended phase (a — 0.5) and one in the localized phase 
[a — 2). We now address in a more general way the 
question of the time-fluctuations of the one-body Green's 
functions'^^ . Since Hhop conserves the total number of 
fermions, the anomalous Green's functions J-mi{t) are 
always zero. Concerning Gmiit), using the expansion 



/" = I^u 'u-np.Cf, we can express it as: 



out) = {^mUi\m) 



fJ.lfJ.2 



i(cMl-«M2)*G'0 



(23) 



where G°^^^ = (^'olcJ^^Cp^ l*o) is the t = Green's func- 
tion of the normal modes. Assuming the model has no 
single-particle energy degeneracy (i.e., if e^j^ = then 
fJ-i = only the diagonal terms with fii — fi2 will 
contribute to the infinite-time average of Gmiit)- 



(24) 



Hence, the time-fluctuations of Grni{t) will be given by: 



SGml{t) 



E 



(25) 



Now we calculate 5g^^ by squaring the previous expres- 
sion and taking the infinite-time average. If we as- 
sume there is no gap degeneracy (i.e., if e^^ — C/^st^mi ~ 
~ £p4#M3 t^"3n ^1 = /i3 and ^2 = t^i) we arrive atPl: 



E 



^1^*2 I 



(26) 



which turns out to be a very useful and sharp tool in the 
analysis of the time-fiuctuations of the Green's functions. 
We will study it in various situations (different quenches 
and different choices of the fermionic operators /„) to un- 
derstand when and why , which is always finite for 
any finite L (due to revivals) , goes to zero in the thermo- 
dynamic limit L — > 00. We will see, in this respect, the 
crucial role played by the weights |wm/ii| l^i^al ^''^'^ by 
localization. 

Let us consider first real-space Green's functions 
Gj-^j^(t), which for ji — j2 = j, have a simple physi- 
cal meaning: the expectation value at time t of the local 
density at the site j. From Eq. (26 1 we get: 



Jl J2 ^ 

Ml#M2 



-'j2/i2 I 



1(^0 I 



(27) 



where is the real-space wave-functions of the eigen- 
state fj,. Depending on a, the eigenfunctions are ei- 
ther localized (for a > 1) or extended (for a < l)p2l. 
To pin-point this, we could monitor the IPR of the fi- 
eigenstate, IPR(/z) = J2j it turns out, all 

eigenstates behave in the same way for the present model: 
either all localized, IPR(/x) > for L — 00, or all ex- 
tended, IPR(/i) ^ 1/L for L — >■ 00, without any mobility 
edge. For that reason, we can just monitor the average 
IPR defined as: 



IPR 



iEipi^(M) = iEE 



(28) 
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shown in the inset of Figure [I (top) for a = 0.5 ("Ex- 
tended" points) and a = 2 ("LocaUzed" points). Corre- 
spondingly, the main panel of Figure |2] (top) shows the 
average value oi S'q , . as a function of the chain size L for 
both choices of a. Notice that when eigenstates are ex- 
tended, the weights \uj^^^ r - 1/L2 in Eq. ^ 
can be essentially taken out of the sum. But, without 
weights, it is a simple matter to show that: 



1 

1? 



E 



I 



< 



L2 



TVf 
L2" 



(29) 



where N-p is the total number of fermions in the initial 
state. Hence, . . is expected to go to zero as 1 /L when 
quenching towards a phase with extended eigenstates 
{a < 1), as indeed found numerically. On the contrary, 
weights are of paramount importance when quenching to 
a phase with localized eigenstates {a > 1), because they 
move the important contributions to Sq, . from the aver- 



m 



age |G°^^2 which is of order 1/L, to rare large valuei 
leading to a finite Sq. . which is rather insensitive to the 
size L. 

Consider now the Green's functions in momentum 
space Gk,k2{t) = (*(^)|cIiCfe2 representing, for 
ki — k2 — k, the expectation value at time t of the 



momentum distribution. Since 



'''^n\l^L, the 



Gkxk-i^ty^ are straightforwardly related to the Gj^j^itys 
through a double summation on ji and j2 with oscil- 
lating phase factors e^'^^^^^~^^^^\ However, for L oo 
these are infinite sums, and this might change the be- 
havior of the time-fluctuations. We now show that, even 
when Gjjj2 {t) has persistent time-fluctuations, the cor- 
responding Gkikiit) averages them out, and Sq^ ^ -> 



for L — > oo. Indeed, using Eq. (26) 



E 



r° I 



(30) 



where 



/V L are the Fourier transforms 



of the real-space wave- functions . Figure p| (bottom) 
shows that Sq^^ — ^ for L — >■ oo, regardless of the value 
of a. Notice that, perhaps counter-intuitively, the aver- 
age IPR in fc-space 



IPRk-spacc = 7 E E \'^kti\'* 



(31) 



always decreases as 1/L for both a < 1 (extended real- 
space wave-functions) and a > 1 (localized real-space 
wave-functions). The extended case a < 1 is particu- 
larly intriguing, because one would expect that an ex- 
tended real-space wave-function should look "localized" 
in momentum space, i.e., composed of a small number 
of fc-waves. This expectation, quite reasonable for ordi- 
nary extended states of non-disordered systems, is not 
correct, in general, in presence of disorder, as quite evi- 
dent from Figure [sjj A, B). A simple example demonstrates 
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Spinless fermions - Real space 




Spinless fermions - Momentum space 




^ (bottom) for 

p as a function 
Here G,- 



FIG. 2: Average value of Sq.. (top) and Sq^ 

the disordered long-range hopping model Hhc 
of the chain size L for different values of a. 
(*(f)!c]c,|*(t)) with j = L/2, and Gkkit) = {*(t)|c+Cfe|*(t)) 
with fc = 0. The data are obtained starting from the ground 
state of a clean fermionic chain with nearest-neighbor hopping 
and quenching to Shop with different values of a. The "Lo- 
calized" points are for a — 2, while the "Extended" points 
for a = 0.5. We used 50 realizations of disorder. In all 
cases we report the usual (arithmetic) mean (the error bar, 
when visible, is the standard deviation), except for 5q.. in 
the localized phase, where we plot the median (the geomet- 
ric mean, see Sec. Ill for details). In the inset, the IPR (see 
Eqs. (281 and (Sli) as a function of size. Notice how for 
all cases ("Localized" and "Extended" quenches) the eigen- 
states of the final Hamiltonian are "extended" in reciprocal 
space and consequently 5q^ ^ goes to zero for L — >■ oo. For a 
smoother size-scaling, each disorder realization of the largest 
L generated is employed, by removing the same amount of 
sites from the two edges, to generate realizations for smaller 
L. 



the crucial role played by disorder. Consider a toy real- 
space extended wave- function with Uj = Wj/\/L where 
Wj = ±1 is a random sign on every site. Without Wj, the 
momentum space function Uk = '^IjUjC^^^^ /\/L would 
be localized, with a peak at fc = 0. When Wj is accounted 
for, becomes extremely irregular but extended over 
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FIG. 3: (A) Plot of luj^ip versus the site-index j, where 
^ IS £L typical extended eigenstate of Huop for a — 0.5 
and L = 2048, with energy in the middle of the band. 
(B) The corresponding momentum space |ufc^|'^, with Ukp, = 



quite clearly extended. (C) Plot of \u]^\'^ for 

a toy extended wave-function Uj = Wj/\/Jj, where Wj = ±1 
is a random sign. 



all k points, see Figure [sjjC), indeed with a strong re- 
semblance to the actual momentum space wave-function 
of Figure SJB). In presence of disorder, therefore, being 
"extended in real-space" does not imply a sharply defined 
momentum. Effectively, therefore, going to momentum 
space averages-out persistent time-fluctuations which are 
seen in real space when eigenstates are localized, an effect 
akin to "self-averaging" of extensive quantities in disor- 
dered systemJSl, 



B. Many-body observables and failure of GGE 



When a > 1 (localized phase) the GGE ability in describ- 
ing time-averages of many-body operators is not guaran- 
teed, because Gj^j^{t) have persistent time-fluctuations 
and Wick's theorem is of no help. Here we will show that, 
when a > 1, GGE fails in predicting the spatial density- 
density correlations pj-^ ~ hj^hj^. To see this, we com- 
pare, see Eqs. (21 1- (22), the time-average oi pj^j^{t) (i.e., 
a time-average of a sum of products of G's) with the 
corresponding GGE average (a sum of products of time- 
averages). In Figure |4] we plot Pj^j^it), together with 

its running-time average (i.e., t^^ Jo ^^Viii2 (^')) ^'^d the 
GGE average for two chain sizes and two values of a. 
First we notice that, as in the case of Gj^j^ (0: increasing 
the size L in the delocalized phase (a < 1) strongly de- 
creases the time-fluctuations, which are, on the contrary, 
unaffected by L in the localized phase (a > 1). The 
second feature emerging from Figure [4] is that, while in 
the delocalized phase the time-average tends to the GGE 
value, there is a marked and clear discrepancy between 
the two in the localized phase. 
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From the general analysis of Sees. IV]|V and the study 
of the one-body Green's functions of Sec. |VI A[ we can 
conclude that many-body operators involving a finite ex- 
pansion in terms of momentum space operators Ck^ hke 
correlations c\_^Ck-i^c\^c\^, will have time-averages which 
coincide with GGE averages, regardless the value of a. 
The same is true in the delocalized phase (a < 1) for 
many-body operators with a finite expansion in real 
space, because the timc-fiuctuations of Gj-^^j^ (t) vanish. 



FIG. 4: Time evolution of pj-i^j^if), with ji = L/2 
and j2 = ji + 1 (solid line), and its running-time average 
Jq dt' pj-^j^i^') (dashed line), for two different values of a 
and two different sizes. The horizontal dash-dotted line is the 
GGE average for Pjija- The data are obtained using a single 
realization, see caption of Figure [l] for details. 

The difference between the GGE average and the time- 
average can be explicitly computed using the same strat- 
egy (and assumptions) of Sec. VI A 
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{Pjlj'l) GGE (PjU2)timc 

(4i"=Ji)gGe(4cJ2)gGE - (4cJ2)gGE + r^^]^ f ( | Gj, (t ) ^ - Gj, ) Gj, )) 





= lim — 

T^oo T 



dt{\SG,,,,{t)\'-SG,,,At)SG,,,M 



(32) 



where we first used the Wick's expansions of {pj-ij2)GGE 
and (^'(i)|/5j^j2 then used the relationships 



and 



(c]^c]^ )gge (since GGE works for one-body averages), 



and finally made use of Eq. (25) and of the no-gap- 



degeneracy assumption. The result is closely reminis- 
cent of Eq. ( 27 1 for 6^, , , except that now the weights 

have two contributions. Using the relation IG?. „ P — 
|G^2;^iP we can finally re-express Aj^j^ as an explicitly 
positive quantity as follows: 



A, 



^ 2 $Z [l^jlPlI l'"j2A'2l +l'"jlA'2l l%2A'll (^ilAll'"ilM2''^j2Ml'"j2M2 '^■'^■)] I^MlA'2l 

^1/^2 

— 9 l'"jlMl'"i2M2 ~ ^il/i2^j2/il I |^°i/i2l — ^ ■ (^'^) 



This explicitly shows that, apart from the trivial cases in 
which all terms of the summation are zero (e.g., ji = j2 = 



J 



when 



Pji 



hj becomes a one-body operator 



for which GGE works and Ajj = 0), the GGE average 
overestimates the time-average of Pj^j2 (i), as exemplified 
in Figure |4] for a single realization. The question is now 
if Ajjj2 goes or not to zero for L ^ oo. In Figure [5] 
we plot the average value of Ajjj2 as a function of L for 
the same set of quenches presented before. [Aj-jj^jav goes 
to zero in the thermodynamic limit in the "Extended" 
phase, and this confirms our general analysis: the real- 
space Green's functions have vanishing fluctuations in 
the extended phase, and GGE works for many-body 
observables with a finite expansion. On the contrary, 
in the "Localized" phase, [Aji^J^^ ^ cxp([log A^^^J^^) 
remains finite even in thermodynamic limit, and since 
[Ajjjjjav > [Ajij2]Sf (tiy Jensen's inequality), this en- 
sure that also [Aj^j^Jav is finite for L oo. Clearly, 
the persistent time-fluctuations of the Green's functions 
lead to time- correlations between the G,-^ (t) appearing 
in the expansion of Pjij2it), see Eq. (21), which reduce 
the time-average with respect to the corresponding "sum 
of products" of time-averages. 



VII. ISING CHAIN IN TRANSVERSE FIELD: 
RESULTS 



Consider now the case of an Ising chain in transverse 



field, described in Sec. Ill The most important difference 
with respect to the case of spinless-fermions with long- 
range hoppings is the presence of particle non-conserving 
(BCS-like) terms, which mix particles and holes. Be cause 
of that, it is convenient to work with Nambu vectord^^^H^ 
which make the algebra very similar to the previous one, 
with very similar results. In particular, we will show that, 
due to disorder and localization of eigenstates, the time- 
fiuctuations of the local magnetization (t| — 2cjcj — 1 
remain finite; however considering for instance the total 
magnetization = ctJ /L which is extended over the 
whole chain, introduces an infinite summation which ef- 
fectively leads to a "self-averaging" of time-fiuctuations. 

Consider, for simplicity, only real-space Green's func- 
tions. In terms of the real-space fermion operators Cj 



At 



define the Nambu vector ^ = ^ci, . . . , cl, Cj^, . . . , 
Similarly, in terms of the fermions 7^ which diago- 
nalize Hxy we define a second Nambu vector F = 

T 



,71, 



,7l 



These two vectors are con- 



nected through the relation = UF where U is a 
2L X 2L unitary matrix which describes the Bogoliubov 
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7^, and with energy > 0) there is a corresponding 
index p. = fi + L G [L+1, 2L] (associated to 7^, and with 
energy — e^); similarly, for every /i e [L + 1, 2L] we define 
fi — ^ — L (z The time evolution of is simple: 



iHti 



10' 



e""*r^e = e '^'^T^ where the energy = when 
1 < /i < L and = — when L+1 < fj. < 2L. The pairs 
Cj-^Cj2 and cj^cj^ entering in the one-body standard and 
anomalous Green's functions can be both obtained from 
the Nambu pairs 'fj^^jj (where now ji and j2 run from 

1 to 2L): ^]^*,, = E^,^,f/;,;.,C^..M.fj.,f^., and the 
associated (Nambu) Green's function closely resembles 
Eq. ([23l): 



FIG. 5: Average value of Aj^j^ = {pjij2)GGE — {pjij2)tiaie, the 
discrepancy between the GGE average and the time-average, 
for the density-density correlation Pj^j^ = njj 



tion of the chain size, for with ji = L/2 and two values of 
j2 — ji ■ The data are obtained using the same quenches of 
Figure [2j see its caption for details. In the localized phase 
(a = 2) we plot the median [Aj^jg, = exp([log Aj^jJ^^) 
because A^^j^ is there roughly log-normal distributed. 



'-^'Sl,. - (34) 



as a func- where = (^olrj^^r^J^'o) . Assuming, again, no 

energy degeneracy, the time-fluctuations of Gj-^j2 (t) read: 



^A'lM2 



(35) 



rotatiorP^I^ performed to diagonalize H 



XY- 



Notice that the index /i in the Nambu vector F^ runs 
from 1 to 2L: for every index /i G (associated to 



In computing the time-averaged squared-fluctuations of 
Qmi(t) we have to take care of gap degeneracies due to 
the particle-hole symmetry of the spectrum e^. Taking 
due care of that, the value of 5^ turns out to be: 

^3132 



nO 

U2U\ 



^1,32 = E E e^^'--'--'-+'-^*c^..c^2.2c/,...^..^;:,M.e: 
= E i^jiMin^i2A'2n^°iM2r 

A'l#A'2 



E ^JlMl^i2Ai2^JlA'2f^j2?ir l^/^lA'2 



(36) 



where the over-line denotes the infinite-time average, and defined similarly to Eq. (28) 



we used the relation F^, 



The first term is due to 



the cases in which /ii ~ vi and ^2 = ^^2, similarly to 
what found in Eq. ( 26 1 . The second term originates from 



particle-hole symmetry (present even when the system is 
disordered) and occurs when fii — and fi2 — v\- As 
we did for the spinless fermionic Hamiltonian, one can 
show that: 



(37) 



iVIoreover, using the fact that J7^[/ = 1 we can repeat the 
same observations presented in Sec. |VI A[ if the eigen- 
states are delocalized (clean chain case) Sg, _ goes to 

zero for L — > (X), while Sj<, remains finite when the 
eigenstates are localized. For any finite disorder ampli- 
tude e, the Hamiltonian Uxy has always localized states. 
Indeed, as shown in the inset of Figure[6] the average IPR 



2L 2L 



IPR 



Isine 



(38) 



is finite for L 00. This localization leads to persistent 
time-fluctuations: 6g, . > in the thermodynamic limit. 



This is shown in Figure |6] where we plot the average value 
of (5^z = 45g . . for a quench from the ground state of a 

clean Ising chain at the critical point (e = 0, 7 = 1, 
Jj = 1, /ij = 1) to a disordered Ising chain at the infinite 
randomness critical point (e = 1, 7 = 1, Jj e [0,2], 
/i,€[0,2]j3sl. 

We now show that, while each c|cj has nonvanish- 
ing time- fluctuations, the average magnetization per site 
rhz — cF^j /L ^2 Y^^=i j / has vanishing time- 
fluctuations for large L, due to cancellations reminiscent 
of "self-averaging" in extensive observablepl'. Indeed, 
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FIG. 6: Value of 5%, with 6 the local transverse magnetiza- 
tion 6"! or the total magnetization rh^, as a function of the 
chain size. The initial state is the ground state of a clean 
Ising chain at the critical point (e = 0, 7 = 1, Jj = 1, hj — 1) 
and the final Hamiltonian is a disordered Ising chain at the 
infinite-randomness critical point (e = 1, 7 = 1, Jj £ [0,2], 
hj € [0,2]). Data obtained with 50 disorder realizations. For 
a smoother size-scaling, each disorder realization of the largest 
L generated is employed, by removing the same amount of 
sites from the two edges, to generate realizations for smaller 

L. For O = cTj we plot the median [(5^^]^/ ~ sxpnlogS^z ) 

J L J J av 

because S'^z is roughly log-normal distributed. 



Smz{t) — 2^ ■ SGjj{t)/L, which implies that: 



6mz{t)= ^ 



i.e., an expression entirely similar to Eq. ( 35 ) for SQj-^j^ (t) 
except for the vireight U*_^^_^ ^^2^2 which is now replaced by 

the averaged weight w^^^^ = '^J2j=i Ujfj,^Ujf_,^/ L. With 
the same steps done to obtain Eq. (36), we finally get: 



I^A'lM2 I 



-E 

Pl7^P2 
Pl7^P2 



* 1/^0 I 
2^7n7i2 I='miM2 I 



which is definitely different from Eq. (36), the site av- 
erage having been performed before taking the squared 
time-fluctuations. In Figure[6]we plot S'^^ , averaged over 
disorder realizations, as a function of the chain size L: we 
clearly see that, even if the eigenstates of the Hxy are lo- 
calized, the time-fluctuations of decay, and ^ 
for large L. This behavior for rhz is similar to that of the 
Green's functions Gkik2{t) for the disordered long-range 
hopping fermions analyzed in Sec. |VI A[ where the infi- 
nite site summations lead to a cancellation of the time- 
fluctuations of the various terms. 



VIII. DISCUSSION AND CONCLUSIONS 

Let us discuss some of the most relevant recent pa- 
pers appeared in the literature, in the light of what 
we have presented in our paper. A detailed analysis 
of the validity of GGE averages for quantum quenches 
where the final Hamiltonian was integrable, disorder-free 
and translationally invariant (the ID quantum Ising/XX 
spin chains and the Luttinger model) has been made by 
Cazalilla et al^, showing that, for a general class of ini- 
tial states [vPo), the time-fluctuations of the one-body 
Green's functions vanish and the GGE averages are cor- 
rect, in the thermodynamic limit, for both local and non- 
local observables. These results are in complete agree- 
ment with what we have shown here, since homogeneous 
Hamiltonians have extended eigenstates and the time- 
fluctuations of the one-body Green's functions decay for 
i — >■ 00. In our study, we have extended the analysis of 
Ref. [3T] to quantum quenches with a final Hamiltonian 
H which is disordered: we have shown that the local- 
ization properties of H are crucial for the relaxation of 
time-fluctuations and, ultimately, also for the validity of 
the GGE. In particular we showed that, for one-body op- 
erators, infinite-time averages are always reproduced by 
GGE, while, for many-body operators, the localization 
of eigenstates of H and the ensuing absence of relaxation 
of one-body real-space Green's functions are in principle 
dangerous for the validity of GGE. 

Quantum quenches with integrable Hamiltonians hav- 
ing a transition between extended and localized states 
have been analyzed in a recent work by Gramsch et al^^. 
They have studied the Aubry-Andre model^^ for hard- 
core bosons in a one-dimensional quasi-periodic poten- 
tial, H = Y.j=l{^k+l + H-c.) AX^^^i cos(27rcrj)nj, 
where S| (bj) creates (annihilates) a hard-core boson at 

site j, fij — bjbj, a is an irrational number and A is 
the strength of the quasi-periodic potential. This model 
can be diagonalized by mapping it, through Jordan- 
WigneJSSl^ onto a non-interacting spinless fermion chain 
with the same potential. The quasi-periodic on-site po- 
tential, in absence of a true disorder, is able to induce a 
transition to a phase with localized eigenstate^^ at a fi- 
nite strength of A = 2. Ref. HOI considered, among others, 
two operators that are particularly relevant for our dis- 
cussion: the local density of bosons hj , a one-body oper- 
ator in terms of Jordan- Wigner fermions, and the boson 

momentum distribution TOfe = X^jus ^''^^"'^~"'^^^1i^J2/-^j 
which is, on the contrary, a many-body operator when 
written in terms of fermions, because of a Jordan- Wigner 
strinsBSl. The numerical results of Ref. 40 show that 



when the eigenstates of the final Hamiltonian are ex- 
tended (A < 2), the time-fluctuations of both hj and m^, 
vanish, and the GGE predicts the time-averages quite 
well, consistently with our analysis; when the eigenstates 
of the final Hamiltonian are localized (A > 2) the sit- 
uation is more complex: the time-fluctuations of hj do 
not relax but GGE predicts well the time-average (again 
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consistently with our analysis of Sec. IV I, while the time- 



fluctuations of TTik appear to vanish, but GGE seems to 
fail. The failure of GGE in predicting time-averages of 
a many-body observable like rhk when persistent time- 
fluctuations of the one-body Green's functions are at play 
(localized phase, A > 2) is perfectly in line with our re- 
sults (see Sec. VI B). What is definitely beyond our anal- 



ysis, but not in contradiction with it, is the fact that the 
time-fluctuations of rhk relax for large t: this is likely 
an effect of cancellation of fluctuations due to the sum- 
mation of many terms, similar to what wc have found 
for extensive operators (see Sec. [VII ) or for momentum 
space Green's functions (see Sec. VIA). 

Another paper quite relevant for our study is that of 
Khatami et al.^^\ where they analyze quenches with a 
final Hamiltonian similar to our Hhop, Eq. (l8|, supple- 
mented by an interaction term V ^^(ni — l/2)(ni+i — 1/2) 
which definitely breaks integrability. Two comments are 
here in order. First, as shown in Ref. [311 interactions 
do not change the picture dramatically: numerically, a 
metal-insulator transition occurs around a ~ 1^1.2, with 
a quite clear metallic phase for a ^ 1, and an insulat- 
ing one for a > 1.2. Second, by comparing after-quench 
time-averages with the micro-canonical average for the 
momentum distribution function c^Ck and the density- 
density structure factor Y^^i e^''^'-~"^^nihm/ L, Ref. IMI 
shows that quenches in the metallic phase (a < 1) are 
well described by the micro-canonical ensemble, while 
thermalization appears to break-down when quenching 
to the insulating phase {a > 1.2). These results are 
definitely in line with what we have found, and suggest 
that, independently of the integrability of the Hamilto- 
nian, the localization properties of H are crucial for the 
after-quench thermalization. 

Finally, let us mention a technical point related to the 
relaxation of time-fluctuations of general many-body op- 
erators of nonintegrable model^SSIMI with the technique 
we have used for one-body Green's functions of free- 
particle models. In principle, one could compute for 

a general operator O, starting from Eq. ([2|. However, 
to make progress, one would need to stipulate something 
about gap degeneracies in the many-body spectrum, i.e., 
Ea — Ep — Ea' — Ep> with a ^ a' and 13 ^ 13' (apart 
from the trivial case a = /3 and a' = f3'). (The as- 
sumption of absence of gap degeneracies is often used in 
the literature, and dates back to the original paper of 
von Neumann^ who, however, carefully stipulates it to 
hold only within each micro-canonical energy shell, and 
not for the many-body spectrum at large.) The condi- 



'Ml 'M2 
En 



Ep> 



and 



En 



tion of absence of gap degeneracieJ^ is clearly unten- 
able for models with non-interacting quasi-particles: you 
can produce an exponentially large number of many-body 
states \a') and \/3') whose spectral gap Eai — Ep' coin- 
cides exactly with Ea — Ep: simply operate on \a) and 
1/3) by applying, in identical fashion, an arbitrary number 
of particles and/or holes, \a') 

W) =lUlU---lUP)- Thenii;^, 

+ + • • • e^^ because quasi-particles do not inter- 
act, and therefore E^' — Epi — Ea — Ep. One might 
argue that this proliferation of exactly degenerate gaps 
is a peculiarity of models with non-interacting quasi- 
particle^^ interaction effects might change the picture 
completely. This is definitely an interesting point, which 
deserves further studies, but certainly also a very hard 
one, because the combination of disorder and interac- 
tions makes the analysis highly nontrivial. Neverthe- 
less, let us mention the following simple argument. Sup- 
pose that quasi-particles interact, but there is still an 
exponentially large (in the number of particles N) num- 
ber of states with very small spectral gap differences A: 
Ea — Ep = Ea' — Ep' + A. Then, all these spectral gap 
quasi-degeneracies will appear, effectively, as true degen- 
eracies until a time T h/ A is reached, and that time 
might indeed be very large. 

In conclusion, we have analyzed after-quench relax- 
ation and thermalization issues for "solvable" models 
with disordered fermions (or spins) in one dimension. 
Several points are still open and deserve further stud- 
ies. Let us just mention two of them: What is the role 
of a mobility edge, separating localized from extended 
states, in the after-quench relaxation? What is the role 
of interactions in combination with disorder in the out- 
of-equilibrium dynamics of closed quantum systems? 
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